Expanding (3+l)-dimensional universe from a Lorentzian matrix model 
for superstring theory in (9+l)-dimensions 
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We reconsider the matrix model formulation of type IIB superstring theory in (9+l)-dimensional 
space-time. Unlike the previous works in which the Wick rotation was used to make the model 
well defined, we regularize the Lorentzian model by introducing infrared cutoffs in both the spatial 
and temporal directions. Monte Carlo studies reveal that the two cutoffs can be removed in the 
large- iV limit and that the theory thus obtained has no parameters other than one scale parameter. 
Moreover, we find that three out of nine spatial directions start to expand at some "critical time" , 
after which the space has SO (3) symmetry instead of SO (9). 
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Introduction. — One of the most fundamental ques- 
tions concerning our universe is why we live in a (3+1)- 
dimensional space-time, and why the universe is expand- 
ing. The aim of this Letter is to provide some evi- 
dence that these facts can be derived from a nonper- 
turbative formulation of superstring theory in (9+1) di- 
mensions based on matrix models. Motivated by re- 
cent developments in understanding the dynamics of the 
Euclideanized model, we study the SO(9, 1) symmetric 
Lorentzian model nonperturbatively without Wick rota- 
tion. Our Monte Carlo results demonstrate, among other 
things, that three out of nine spatial directions start to 
expand in the early universe. We expect that what we 
are doing here is essentially a first-principle calculation of 
the unified theory including quantum gravity. This may 
be contrasted with the quantum cosmology in the early 
1980s that aimed at describing the birth of the universe 
l] within the mini-superspace approximation. More re- 
cently, a nonperturbative approach to quantum gravity 
has been pursued using the causal dynamical triangula- 
tion [2j. For earlier works that put forward the idea to 
use matrices for cosmology, see Refs. 0,0]. 

Matrix model for super strings. — Superstring theory 
not only provides a most natural candidate for a consis- 
tent theory of quantum gravity but also enables a unified 
description of all the interactions and the matters. A cru- 
cial problem is that we do not yet have a well-established 
nonperturbative formulation, which would be needed in 
addressing dynamical issues such as the determination of 
space-time dimensionality 

In the 1990s, there was remarkable progress in under- 
standing the nonperturbative aspects of superstring the- 
ory based on D-branes. Most importantly, it was noticed 
that large-A^ matrices are the appropriate microscopic 
degrees of freedom which are useful in formulating su- 
perstring theory in a nonperturbative manner [7H9J. In 



particular, the type IIB matrix model was proposed as a 
nonperturbative formulation of type IIB superstring the- 
ory in ten-dimensional space-time [8J . It was also realized 
that the five types of superstring theory in ten dimen- 
sions are just different descriptions of the same theory. 
Therefore, it was speculated that the type IIB matrix 
model actually describes the unique underlying theory, 
although it takes the form that has explicit connection 
to perturbative type IIB superstring theory. 

In the type IIB matrix model, the space-time is rep- 
resented dynamically by the eigenvalue distribution of 
ten bosonic N x N traceless Hermitian matrices [Icj . 
So far, the dynamical generation of four-dimensional 
space-time has been discussed exclusively in the Eu- 
clideanized model. Indeed the spontaneous symmetry 
breaking (SSB) of SO(10) down to SOU) was suggested 
by the Gaussian expansion method [HI [l2j ]. Recently, 
systematic calculation of the free energy has been per- 
formed for SO(d) symmetric vacua with 2 < d < 7, and 
it is found that d = 3 gives the minimum (l3 |. Further- 
more, the ratio of the space-time extent in the extended 
directions to that in the shrunken directions is shown to 
be finite. These results, if true, suggest the necessity 
for reconsidering the formulation in order to make any 
connection to the real world. 

Matrix model with SO (9,1) symmetry. — Our start- 
ing point is the action S = Sh + Sf , where [8j 



S h = -^tv{[A^,A u ][A^A"- 
St = —ix^fciCT^aplA^p]) 
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with Af, (/i = 0, ••• ,9) and * Q (a = 1, • • • ,16) be- 
ing N x N traceless Hermitian matrices. The Lorentz 
indices [i and v are contracted using the metric r\ — 
diag(— 1, 1, • • • ,1). The 16 x 16 matrices are ten- 
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dimensional gamma matrices after the Weyl projection, 
and the unitary matrix C is the charge conjugation ma- 
trix. The action has manifest SO(9,l) symmetry, where 
and transform as a vector and a Majorana-Weyl 
spinor, respectively. The Euclidean model, which has 
SO(10) symmetry, can be obtained from this action by 
the Wick rotation Aq — iAiQ. A crucial difference is that 
the bosonic part of the action in the Euclidean model is 
positive definite, whereas in the Lorentzian model it is 



tr {F liV F ttv ) = -2 tr (F Ql ) 2 + tr (F^f 



(2) 



where F^ v = -i^,^] are Hermitian matrices, and 
hence the two terms in ^ have opposite signs flij ]. 

We study, for the first time, the Lorentzian model non- 
perturbatively based on the partition function 



Z = J dAd-V e iS = J dAe lSb PfM(A) 



(3) 



where the Pfaffian Pf M. (A) appears from integrating out 
the fermionic matrices Note that in the Euclidean 
model, the Pfaffian is complex in general, and its phase 
plays a crucial role in the aforementioned SSB of SO(10) 
symmetry 16|, [17J . On the other hand, the Pfaffian in the 
Lorentzian model is real. Therefore, the mechanism of 
SSB that was identified in the Euclidean model is absent 
in the Lorentzian model. 

In the definition ([3]) , we have replaced the "Boltzmann 
weight" e~ s used in the Euclidean model by e %s . This is 
theoretically motivated from the connection to the world- 
sheet theory [8|. The partition function ([3]) can also 
be obtained formally from pure M = 1 supcrsymmetric 
Yang-Mills theory in (9 + 1) dimensions by dimensional 
reduction. Note, however, that the expression ([3]) is ill 
defined and requires appropriate regularization in order 
to make any sense out of it. This is in striking contrast 
to the Euclidean model, in which the partition function 
is shown to be finite without any regularization 1|| 19| . 

It turns out that the integration over Aq is divergent, 
even if we fix -itr (Ai) 2 to a constant. In order to cure 
this divergence, we introduce a constraint 



^tr(A ) 2 < K ltr(^) 2 , 



(4) 



which is invariant under the scale transformation A^ 
pA^. Note that this constraint generically breaks SO(9,l) 
symmetry down to SO(9). However, it turns out to be 
equivalent to imposing ^ after "gauge fixing" the boost 
symmetry by requiring that -^-tr {Aq) 2 with A^ ~ 0^ U A U 
be minimized with respect to O € SO(9, 1). In this sense, 
the constraint actually respects the SO(9,l) symmetry. 

Let us note that e lSb in the partition function <j3j> is 
a phase factor just as in the path-integral formulation of 
quantum field theories in Minkowski space. As is com- 
monly done in integrating oscillating functions, we intro- 
duce the convergence factor e -6 ' 5 * 5 ' and take the e — > 
limit after the integration. 



The partition function can then be rewritten as 
Z = J dA J o dr $ (jj tT ( A *) 2 - r ) e tSb - £|Sb| PLM , 

where the integration over A^ is assumed to be restricted 
by the constraint Q . After rescaling the variables A^ — » 
r 1 / 2 A^ in the integrand, we integrate over r and get 



(5) 



which diverges for 5b = 0. In order to cure this diver- 
gence, we introduce a constraint 
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tr (AA 2 < L 2 



(6) 



before the rescaling. Then the integration domain for r 
becomes [0,L 2 ], and (O is replaced by /(5b), where f(x) 
is a function with a sharp peak at x = 0. Thus we arrive 
at the model 

Z=JdAf (±tr (F^F^) - Cj PiM(A) 

x^ltr^-l^^-itrCAo) 2 ) , (7) 

where 8(x) is the Heaviside step function. The constant 
C should be set to zero according to our derivation. If we 
consider the C < case, the model {?]) may be viewed as 
the matrix model motivated from the space-time uncer- 
tainty principle [20j | with the regularization in the second 
line, which we find to be necessary. Since the Pfaffian 
Pf JVl (A) is real in the present Lorentzian case, the model 
([7|) can be studied by Monte Carlo simulation without the 
sign problem. Note that this is usually not the case for 
quantum field theories in Minkowski space. 

Monte Carlo results. — We perform Monte Carlo sim- 
ulation of the model ([7]) with C — by using the Rational 



Hybrid Monte Carlo algorithm [2l|, which is quite stan- 
dard in recent simulations of quantum chromodynamics 
including the effects of dynamical quarks. 

In order to see the time evolution, we diagonalize Aq, 
and define the eigenvectors \t a ) corresponding to the 
eigenvalues t a of Aq (a — 1, • • • , N) with the specific 
order t\ < ■ ■ ■ < tjv. The spatial matrix in this basis 
(t a \Ai\tb) is not diagonal, but it turns out that the off- 
diagonal elements decrease rapidly as one goes away from 
a diagonal element. This motivates us to define n x n 
matrices A^ b \t) = (<„+ a |A;|£i/+h) with 1 < a, b < n and 
t = i J2a=i tv+a for v = 0, • • • , (N — n). These matri- 
ces represent the 9d space structure at fixed time t [22j |. 
The block size n should be large enough to include non- 
negligible off-diagonal elements. In Fig.[T]we plot the ex- 
tent of space R{t) 2 = ±tr A z (t) 2 for N = 16 and n = 4. 
Since the result is symmetric under the time reflection 
t — > — t as a consequence of the symmetry Aq —Aq, 
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FIG. 1: The extent of space R(t) 2 with N = 16 and n = 4 
is plotted as a function of t for five values of ft. The peak at 
t = starts to grow at some critical ft. 

we only show the results for t < 0. There is a critical ft, 
beyond which the peak at t = starts to grow. 

Next we study the spontaneous breaking of the SO (9) 
symmetry. As an order parameter, we define the 9x9 
(positive definite) real symmetric tensor 

i; i (t) = itr{A i (t)A i (t)} , (8) 

which is an analog of the moment of inertia tensor. The 
nine eigenvalues of Tij(t) are plotted against t in Fig. [2] 
for k = 4.0. We find that three largest eigenvalues of 
Tij(t) start to grow at the critical time t c , which suggests 
that the SO (9) symmetry is spontaneously broken down 
to SO(3) after t c . Note that R(t) 2 is given by the sum of 
nine eigenvalues of Tij (t) . 

Mechanism of the SSB.~ The SSB of SO(9) looks 
mysterious at first sight, but we can actually understand 
the mechanism quite intuitively. Let us consider the case 
in which ft is large. Then the first term of ([2]) becomes a 
large negative value, and therefore the second term has 
to become large in order to make ((2|) close to zero as 
required in ((7]). Because of the constraint -^-tr (Ai) 2 = f , 
however, it is more efficient to maximize the second term 
of ([2]) at some fixed time. The system actually chooses 
the middle point t = 0, where the suppression on Ai from 
the first term of ((2]) becomes the least. This explains why 
the peak of R(t) at t = grows as we increase k. 

Let us then consider a simplified question: what is the 
configuration of Ai which gives the maximum itr {Fij) 2 
with fixed jrtv {Ai) 2 = 1. Using the Lagrange multiplier 
A, we maximize the function G = ti (Fij) 2 — Atr (Ai) 2 . 
Taking the derivative with respect to Ai, we obtain 
2 [Aj, [Aj,Ai]] — XAi = 0. This equation can be solved if 
Ai — xLi for i < d, and Ai — for d < i < 9, where Li 
are the representation matrices of a compact semi-simple 



FIG. 2: The nine eigenvalues of Ty(t) with N = 16 and 
n — 4 are plotted as a function of t for ft = 4.0. After the 
critical time t c , three eigenvalues become larger, suggesting 
that the SO(9) symmetry is spontaneously broken down to 
SO(3). 



Lie algebra with d generators. Clearly d should be less 
than or equal to 9. It turns out that the maximum of 
jjtr (F^) 2 is achieved for the SU(2) algebra, which has 
d = 3, with Li being the direct sum of the spin-^ repre- 
sentation and (N — 2) copies of the trivial representation. 
This implies the SSB of SO (9) down to SO (3). The SSB 
can thus be understood as a classical effect in the k — > oo 
limit. When we tune ft with increasing N as described 
below, quantum effects become important. We have con- 
firmed [23] that the n x n matrix Q = Ai(t) 2 has 
quite a continuous eigenvalue distribution, which implies 
that the space is not like a two-dimensional sphere as one 
might suspect from the classical picture. 

Removing the cutoffs. — It turned out that one can 
remove the infrared cutoffs ft and L in the large- TV limit 
in such a way that R(t) scales. This can be done in two 
steps, (i) First we send k to oo with N as ft = (3 N p 
(jp ~ |) (23[. The scaling curve of R(t) one obtains in 
this way depends on (3. (ii) Next we send f3 to oo with 
L. The two limits correspond to the continuum limit 
and the infinite volume limit, respectively, in quantum 
field theory. Thus the two constraints (Ql, © can be 
removed in the large-TV limit, and the resulting theory 
has no parameter other than one scale parameter. 

Let us discuss the second limit (ii) in more detail. 
We find that the inequality © is actually saturated for 
the dominant configurations. Therefore, one only has to 
make the rescaling A^ H> LA^ in order to translate the 
configurations in the model as those in the original 
partition function. It turns out that R(t) for the rescaled 
configurations scales in (3 by tuning L and shifting t ap- 
propriately. In order to see this, it is convenient to choose 
L so that R(t) at the critical time t = t c becomes unity, 
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and to shift t so that the critical time comes to the origin. 
Then R(t) with increasing /? extends in t in such a way 
that the results at smaller \t\ scale. This is demonstrated 
in Fig. [31 where we find a reasonable scaling behavior for 
N = 16 with k = 2.0, 4.0, 8.0. In fact, supersymmetry of 
the model plays an important role here [23j. 
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FIG. 3: Our data for R(t) 2 shown in Fig. Q] with k larger 
than the critical value are replotted against the shifted time 
t — t c in units of the size of the universe R(t c ) at the critical 
time. 



Summary. — In this Letter we have studied the non- 
perturbative dynamics of the Lorentzian matrix model 
for type IIB superstring theory in ten dimensions. In 
order to make the model well defined, we introduce the 
infrared cutoffs on both the spatial and temporal direc- 
tions. We find that the two cutoffs can be removed in the 
large-TV limit. Moreover, the theory thus obtained has 
no parameters other than one scale parameter, which is 
a property expected for nonperturbative superstring the- 
ory. The SO(9) symmetry breaks down to SO(3) at some 
critical time, and the size of the three-dimensional space 
increases with time. The cosmological singularity is nat- 
urally avoided due to noncommutativity. 

There are a lot of questions that should be addressed in 
our model. One of the most urgent questions is whether 
a local field theory on a commutative space-time appears 
at the low energy scale. A possible way is to calculate 
correlation functions of the Wilson loop operators 24J . If 
this question is answered in the affirmative, we consider 
it very likely that our model really describes the birth of 
our universe from first principles. The next step would 
be to show that the four fundamental interactions and 
the matter fields appear in our universe at a later time. 
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